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Abstract 

In this report, we investigate the synchronization of temporal activity in an electrically 
coupled neural network model. The electrical coupling is established by homotypic static 
gap-junctions (Connexin 43). Two distinct network topologies, namely: sparse random, 
network, (SRN) and fully connected network, (FCN) are used to establish the connec- 
tivity. The strength of connectivity in the FCN is governed by the mean gap junctional 
conductance (//). In the case of the SRN, the overall strength of connectivity is governed 
by the density of connections (5) and the connection strength between two neurons (So). 
The synchronization of the network with increasing gap junctional strength and varying 
population sizes is investigated. It was observed that the network abruptly makes a tran- 
sition from a weakly synchronized to a well synchronized regime when (5) or ([/,) exceeds a 
critical value. It was also observed that the (8, /i) values used to achieve synchronization 
decreases with increasing network size. 
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1 Introduction 



Computational models of neural networks have been found to be useful in characterizing and 
validating hypotheses about how information processing occurs in real nervous systems. For 
example, a pulse coupled neural network (PCNN) model in jTj is capable of replicating tem- 
poral neural activity such as spindle waves, sleep oscillations and sustained spike synchrony. 
However, an important issue that affects the study of network dynamics is the choice of neural 
coupling, pp. Neural coupling is accomplished by synapses which can be broadly classified in 
to (a) chemical synapses and (b) electrical synapses. While several studies consider the former 
type of coupling to be the preponderant way of intercellular communication, recent research has 
provided increasing molecular and functional evidence of the latter. More importantly, |2j and 
3J suggest that neurons could also use electrical synapses to achieve intercellular communica- 
tion. There have also been reports that emphasize the importance of electrical synapses in the 
temporal coordination of neuronal activity, 0j, the generation of high frequency oscillations, |5] 
and the generation of oscillatory activity, 0. 

The primary focus of this brief communication is to quantify the extent of synchronization in 
two general network models of electrically coupled neurons. The electrical coupling is achieved 
with the help of static homotypic gap-junctions (connexin-43). The choice of connexin 43 was 
based on a study [7j that presented molecular evidence for its presence in electrical connections 
between pairs of neurons in the visual cortex and hippocampal regions of the juvenile rat brain. 
However, the methods discussed are generic and can be extended to other static and dynamic 
gap-junctions. Two distinct network topologies namely fully connected network (FCN) and 
sparse random network (SRN) are used to establish the connectivity between the neurons. For 
example, Fig.l illustrates a typical gap-junctional connection between a pair of neurons. In 
FCN, the gap-junctions are assumed to be exponentially distributed with mean conductance 
(//). The exponential distribution was chosen as a possible means to capture the non-uniform 
distribution of gap-junctions in neuronal populations. In the case of FCN, the gap-junction 
strength between every pair of neurons is non-zero. An alternate approach to accomplish the 
distribution of gap-junctions is to assume a sparse-random network with a specified density 
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of connections (5). Unlike the FCN, the coupling strength between any two pairs in SRN is 
either zero or one. SRN can be considered as a special case of FCN, where the gap-junctional 
conductance between a pair of neurons greater than a specified threshold is set to one and those 
lesser than the threshold are set to zero. Thus it might not be surprising to view SRN as a 
quantized version of the FCN. 

While neurons are capable of exhibiting a rich set of firing patterns, we consider a population 
of bursting neurons in this study. Bursting behavior in neurons is considered important because 
bursts increase the reliability of synaptic transmission and provide a mechanism for selective 
communication between neurons, jH] • Planar bursters can be classified based on the bifurcation 
mechanism that leads to the corresponding burst activity (see [§] for a summary). In this study, 
we restrict all the bursters in the population to be of the "square-wave" or "fold-homoclinic" 
type. Burst synchronization in general, consists of two components (a) synchronization of 
spikes within a burst, and (b) synchronization between bursts, 18 . In order to minimize the 
contribution of the former, we studied the envelope v^ nv , i = 1 . . . N of the bursts obtained by 
filtering the membrane potentials Vi,i = 1 . . . N. A representative burst and its corresponding 
envelope which approximates the duration of the burst is shown in Fig. 2. 

In the present study, we show that increasing mean conductance (//) and density of connections 
(5) in the FCN and SRN result in increased synchronization as reflected by the synchroniza- 
tion index (M) (Sec. 3). The synchronization index (M) estimated on the original and their 
envelopes is discussed in (Sec. 4). It is also shown that magnitude of the gap-junction strength 
(8, n) to achieve increased synchronization in FCN and SRN decreases with increasing popula- 
tion size. 

2 The Model 

In this study, we represent a network by three attributes, namely (a) individual neurons, (b) 
the connection between neurons and (c) the pattern of connectivity. 
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2.1 Individual Neuron : Model 

While several models are available to represent a single neuron, we choose the recently proposed 
model by Izhikevich pQ, ^0], jTTj which is described by a set of two coupled ordinary differential 
equations 



In Eqn.fjnj), the variable v represents the membrane potential of the neuron while u represents 
a recovery variable which accounts for the activation of K + ionic currents and inactivation of 
Na + ionic currents, [TO]. The strength of the model is that it can exhibit firing patterns of 
every known type of cortical neuron for various choices of the parameters a,b,c and d, [TXUj . 
Moreover, the model is computationally superior to several other neuronal models, while being 
biologically plausible which makes it attractive for conducting large scale simulations and hence 
its choice in the present study. For a complete summary of the neuro-computational properties 
of this and other neuronal models, refer to [HI] . 

2.2 Connections between neurons: Gap Junctions 

Electrical synapses between neurons unlike chemical synapses, are fast and play a crucial role 
in the synchronization of neuronal activity. A gap-junction link consists of adjacent hemi- 
channels called connexons from neighboring cells. Each connexon being composed of proteins 
called connexins, whose conformational structure dictates their conductance. These junctions 
can be broadly classified into homotypic and heterotypic junctions. The variation of the junc- 
tional conductance with the transjunctional voltage is symmetric for homotypic junctions which 
is attributed to identical connexins. However, an asymmetric variation is characteristic of het- 
erotypic junctions. In this study, we implicitly assume the junctions to be represented by the 



v = OMv 2 + 5v + 140 - u + / 



(1) 
(2) 



u = a(bv 



u) 



with auxiliary after-spike resetting given by 




(3) 
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homotypic Connexin 43 (Cx43-Cx43) which has been suggested to mediate neural communi- 
cation, |7j. However, the methods to be discussed are generic and can be extended to other 
types of gap-junctions. Several models have been proposed in the past to model the dynamics 
of gap-junctions [12], [IB], [E]. In the present study, we choose the contingent gating model 
[T4] . where the conductance of a channel between adjacent cells is given by, 



In the above expression, g res represents the normalized residual conductance of the gap-junction, 
9max the normalized maximum conductance, P Q the open probability of the gap-junction; ao, bo 
are the voltage sensitivity coefficients; V\, v% represent the voltages for half-maximal conduc- 
tance. In Eqn.(|5|). v represents the transjunctional voltage which is the difference between the 
membrane potentials of the two cells connected by the gap-junction. Considering a homotypic 
channel (Cx43-Cx43)with identical gates (as suggested in [7j), the coefficients in Eqn.© have 
to satisfy ao = &o and v\ = t>2, |14j . In this case, the variation of junctional conductance is 
symmetric with the transjunctional voltage which means g n (v) = g n {~v)- Further, the gap- 
junctional conductances are assumed to be static. This implicitly assumes the duration of the 
action potential to be much smaller than the changes in the gap-junctional conductance. The 
low-dimensional nature of the static contingent gating model makes it ideal for the present 
study, hence its choice. The values of the parameters in Eqn.(JSJ) for the choice of homotypic 
connexon pairing (Cx43-Cx43) are presented in Table 1, J3]. With the present model, a dif- 
ferent connexon pairing can be obtained by choosing the appropriate values of the parameters 
{dres, 9max, a o,Vi, bo an d v 2 ) for that pairing from Table 1, [Hj. While several factors affect the 
strength of gap-j unctions, in the present study we represent the cumulative effect by a single 
parameter (S). Thus, the conductance of a gap- junctional channel between two neurons (i, j) 
in the network is given by 



A schematic representation of a pair of neurons connected by a symmetric homotypic gap 
junction is shown in Fig. 1. 



P Q = {1 + e «o(-«-«i) + e bo(v-v 2 )y- 




(4) 
(5) 




(6) 
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2.3 Pattern of Connectivity : SRN and FCN 

In an ideal network, the neurons should be connected to accurately reflect the pattern of neural 
interconnections in the real brain. However, a lack of precise knowledge of the anatomical 
connectivity makes such a construction difficult in general. Here, we consider two models for the 
pattern of interconnections namely (i) a sparse random network (SRN) with fixed conductances 
through out the network and (ii) a fully connected network (FCN) where the conductances 
are exponentially distributed about a mean conductance \x. In the sparse random network 
model, we consider the gap-junctional strength (S) to be fixed between any two neurons, i.e. 
Sij = Sji = So if the pair of neurons are connected. The random network is created by 
generating a symmetric matrix W (N x N) with density 5 such that W has (N 2 5) non-zero 
entries. Here, the density 5 of the SRN is defined as the of connections) /N 2 . Then, a 
connection between a pair of neurons (i, j) is assigned only if W(i,j) ^ 0. If W(i,j) = 0, 
then no connection exists between neuron i and neuron j. Thus, the sparsity of the network 
is controlled by the parameter 8. In FCN, we assume a fully connected structure, that is, 
every pair of neurons in the network is connected. However, the connection strengths are 
assumed to be exponentially distributed about a mean strength /i. In both networks, note 
that the connections are made symmetric to enforce the bi-directional nature of gap-junctional 
couplings. 

2.4 Overall Network Model 

The overall network model is obtained by combining Eqns.^EJ) for the neurons, Eqns. (j4l5|) for 
the gap junctional conductances and summing the contributions of synaptic current injections 
resulting from the appropriate network topology. In complex physiological systems such as 
neural assemblies, discrepancies in ionic exchanges between the neurons and their environment 
render the quantitative dynamics of individual neurons different from one another. Therefore, 
the neural population is modeled as set of N bursting neurons whose intrinsic parameters Cj, di 
and the thalamic input Ii are perturbed so that each neuron has a different burst duration and 
frequency. This is done by setting = 0.02, hi = 0.2 and (q, di) = (—65, 8) + (15, — 6)rj where 
Ti is a random variable uniformly distributed in [0,1]. Thus, the dynamics of the i th neuron 
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in the network is described by 



Vi = 0.04^ 2 + hvi + 140 - Ui + h + 9ij{vj - Vi) 



(7) 



3 

Ui = ai(biVi - Ui) 



(8) 



with auxiliary spike- resetting Eqn.fjHJ) 

3 Measures of Synchronization 

Several measures have been proposed in the past jTH] to determine the extent of synchronization 
in neural systems. These can be broadly categorized into linear and non-linear measures. While 
nonlinear measures may be more appropriate in the present context, their estimation can be 
quite challenging. In this report, we use pair-wise linear correlation (0 < p < 1) and Morgera's 
covariance complexity (0 < C < 1) jTH] to quantify the extent of synchronization between the 
neurons. 

3.1 Pairwise Correlations 

The index {pij) between two neurons (i, j) is given by the linear correlation between their mem- 
brane voltages Vi, Vj respectively. For a system of N neurons, we obtain P = N C2 = N(N—l)/2, 
pair-wise dependencies. Statistically significant pair-wise dependencies were chosen as those 
whose p- values are lesser than a specified level of significance (a = 0.05). The null hypothesis 
addressed is that there is no significant correlation between a given pair of neurons. It should be 
noted in the present study we have P pair-wise dependencies. In statistical literature, the choice 
of multiple testing correction such as Bonferroni correction is often recommended to minimize 
the false-positive rate. The adjusted significance is given by a* = a/ P. For example, in a ten 
element (N = 10) network, we have P = 45 and therefore, a* = a/P = 0.05/45 = 0.0011. In 
the present study, the proportion of significant pair-wise correlations was determined as those 
whose p-values were lesser than a divided by P. 
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To determine whether the distribution of the (P) pair-wise dependencies across any two states 
of activity for the network is statistically significant (i.e. a = 0.05), we used parametric fit- 
test) and non-parametric (Wilcoxon-ranksum) tests. Unlike the parametric test (t-test), non- 
parametric tests (Wilcoxon-ranksum) do not assume a normal distribution of the values and 
hence unbiased. However, non-parametric tests are based on ranks and not the original values, 
which is a limitation. Thus using a combination of parametric and non-parametric tests can 
minimize spurious conclusions that are an outcome of a particular test's assumption. 

3.2 Morgera's Covariance Complexity 

Morgera's covariance complexity is a linear measure and has been used recently for the analysis 
of brain signals JH], Here, we use this measure to define a synchronization index as follows. 
After simulating the activity of the network for a given interval (T) using a time step of h = 0.5 
ms, the resulting M = T/h discrete time points of the membrane voltages of the N neurons in 
the system is represented as a matrix T MxN where (N << M). Singular value decomposition 
(SVD) of T yields N eigenvalues which explain the variance along the orthogonal directions in 
N-dimensional space. The normalized variance along the i th component is given by 



i=l A 

Then, Morgera's covariance complexity C is given by. 



^* = ^dhf * = 1 --- N ( 9 ) 



1 p=N 

C = -^Y,^ g a p (10) 

The synchronization index (M) is given by M = (1 — C) and lies in the closed interval [0,1]. A 
minimum value of (M = 0) is obtained in the case of random behavior, whereas a maximum 
value (M = 1) is obtained in the case of perfect synchronization. However, due to external 
noise and nonlinear effects, the estimate of M deviates from these extreme values. 
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4 Simulation Results 



The effect of increasing gap-junctional conductance on the synchronization was determined for 
a given population of neurons. Two distinct network topologies namely, FCN and SRN were 
considered. The synchronization is studied with increasing /i in case of FCN and increasing 
5 in case of the SRN. The synchronization index (M) was estimated on the actual membrane 
potentials (vi,i = 1 . . . N) and their envelopes (t> 4 em ', i = 1 . . . N). The choice of the envelopes 
was in order to minimize the effect of spikes within bursts on the synchronization index. In 
order to determine the effect of varying population size on the synchronization, we investigated 
three different populations namely (N = 10,50 and 100). 

The estimate of M determined on 20 independent trials with varying gap-junctional conduc- 
tance and population sizes is shown in Figure 3. In the case of FCN, the synchronization index 
saturates around (M ~ 0.9) for parameter (// = 40, 20 and 5) and population sizes (N = 10, 50 
and 100), Figure 3 (top). An abrupt transition in the synchronization index is observed at a 
critical density of 5 ~ 0.2 for all three network sizes, Fig.3 (bottom). To generate a quantized 
version of the FCN, we constructed SRN, where the connection strengths So was fixed at 40, 
20 and 5 for populations (N = 10, 50 and 100) respectively. While the connection strength (S ) 
is fixed between connected pairs of neurons, their density 5 (Sec 2.2, Eqn.(6)) which controls 
the number of connections is gradually increased. It should be noted that while 5 G [0, 1] is a 
normalized measure of network density, the actual number of connections (N 2 5) varies in the 
range [0, iV 2 ] for a network of size N. For clarity, we show the variation of the synchronization 
index (M) with respect to (S x S — S x ^number of connections/ N 2 ) which we shall denote 
by S. Note that S quantifies the average connection strength of the SRN. As with the case of 
FCN, we find the network synchronizes with increasing S. The synchronization is achieved at a 
magnitude similar to FCN. As expected, increasing (//) or S leads to increasing synchronization 
of the network. Interestingly, the magnitude of the parameter (//, S) required to achieve syn- 
chronization decreases with increasing population size in the FCN and SRN respectively, Fig. 3. 

While the estimate of M obtained on the membrane potentials and their envelopes converge 
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to a similar value with increasing (S or /i), the abrupt transition is better reflected by those 
estimated on the envelopes, Figure 3. This behavior was consistent varying population sizes 
(N = 10,50 and 100). 

Parametric (t-test) and non-parametric (Wilcoxon-ranksum) statistical tests were used to de- 
termine statistically significant change (a = 0.05) in the distribution of the pair-wise linear 
correlation before and after synchronization. The parameters (S = 0) and (fi = 0) were used 
to represent the networks SRN and FCN before synchronization. The (/i) and (S) values 
after synchronization for the various population sizes (N = 10, 50 and 100) were chosen as 
(/!, S = 40, 20 and 5) from Figure 3. The distribution of the pair-wise linear correlations before 
and after synchronization for the various population sizes for FCN and SRN is shown in Figure 
4. The distributions were statistically significant after Bonferroni correction. This result was 
consistent for all three population sizes. 

5 Conclusions 

In this brief report, we studied the synchronization of temporal activity in a neural network 
with electrical synaptic couplings. The neuronal parameters were randomized to yield different 
burst durations and frequencies but constrained to restrict their behavior to square-wave type 
bursting. The electrical couplings were established by a generic model to represent the homo- 
typic, bi-directional and static gap-junctional couplings formed by Connexin-43. Two network 
topologies namely, (a) sparse random network (SRN) and (b) a fully connected network (FCN) 
were used to model the pattern of connectivity. The overall mean conductance of the network 
was governed by S the density of interconnections and So the fixed connection strength in SRN, 
and \i the mean conductance in FCN. Synchronization of burst activity was quantified using (a) 
the distribution of pairwise linear correlations and (b) Morgera's covariance complexity. Nu- 
merical simulations were conducted to study the synchronization of the network with increasing 
coupling strength and varying population sizes. 

We observe that in both networks, an increase in the mean gap junctional coupling strength 
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(S, fj) results in an increased degree of synchronization. In addition, we observe that the net- 
work exhibits an abrupt transition from a poor to highly coordinated regime at a critical value 
of these parameters. Further, while the magnitude of the parameter \x required to achieve 
synchronization decreases with increasing population size in FCN, the abrupt transition to in- 
creased synchronization seems to occur at a network density of 5 ~ 0.2 for all three network 
sizes in case of the SRN. 

Our results in general, support the prevailing notion that gap-junctions promote synchrony in 
large neural assemblies. The computational studies show that in addition to coupling strength 
of the gap-junctions, the density of their connections in the two topologies (SRN and FCN) also 
plays an important role in promoting synchrony. Moreover, the results suggest the existence of 
critical parameters at which a large neural network could tune, or transition to a high degree 
of synchronization. 
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Figure 1: A pair of neurons coupled by a symmetric homotypic gap junction such that gij(vij) = 
gji(vji) where v = Vij = —Vji is the transjunctional voltage. The tonic activation currents are 
Ii and Ij for neurons i and j respectively. 
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Figure 2: A representative burst (solid line) and its corresponding 'envelope' (dotted line). The 
envelope approximates the duration of the burst. 
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Figure 3: Variation of the synchronization index (M) with parameters ji and 5 for the fully 
connected network (top) and sparse random network (bottom). The vertical lines represent the 
variance about the mean value (dots) for twenty independent trials. Estimates of (M) obtained 
on the original waveforms Vi,i = 1 . . . N (solid lines) and the envelopes vf nv , i = 1 . . . N (dotted 
lines) for three different population sizes N = 10, 50 and 100 is also shown. 
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Figure 4: Distribution of the pair-wise linear correlations estimated before and after synchro- 
nization for FCN (top) and SRN (bottom). Pair- wise correlation before synchronization for 
FCN and SRN were obtained by setting (/x = 0) and (S = 0), represented by hollow bars. 
Distribution of the pair-wise correlation after synchronization for parameters (/x = 40, 20 and 
5, top) and (S = 40, 20, 5, bottom) for population sizes (N = 10, 50, 100) is shown in (a, b, c) 
and (d, e, /) respectively 
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